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Abstract 

Two different theoretical formulations of the finite temperature effects have been 
recently proposed for integrable field theories. In order to decide which of them 
is the correct one, we perform for a particular model an explicit check of their 
predictions for the one-point function of the trace of the stress-energy tensor, a 
quantity which can be independently determined by the Thermodynamical Bethe 
Ansatz. 



1 Introduction 



Finite temperature correlation functions are important quantities for many applications 
of both theoretical and experimental interest (see, for instance A special class of 
quantum field theories is provided by the two-dimensional integrable models, which can 
be exactly solved by means of bootstrap methods 0, ^ ^, ^. For these models, two 
different formulations of finite temperature effects have been recently discussed in the 
literature: the first is due to LeClair and Mussardo 0, the second has been proposed by 
Delfino [0]. Although the two formalisms coincide if applied to the trivial cases of free 
quantum field theories, however they drastically differ once used to deal with interacting 
theories. To determine which of the two is the correct one we decide to compare their 
predictions versus a quantity which can be independently determined. This is the case of 
the finite temperature one-point function of the trace of the stress-energy tensor which can 
be computed by the Thermodynamical Bethe Ansatz (TBA) f^. As we will show below, 
the proposal by LeClair and Mussardo exactly matches the low-temperature expansion of 
this quantity whereas the proposal by Delfino fails at order (9(e~^"*''). Before presenting 
the explicit calculations, let us briefly discuss the main features of the two different finite 
temperature formalisms. 

2 LeClair— Mussardo Formalism 

This formalism, discussed in Q, combines together physical principles coming from two 
different areas: the Thermodynamical Bethe Ansatz and the Form Factor Approach. It 
originates from an interpretation of the expression of the free energy - as determined 
by the TBA -, in terms of quasi-particle excitations with respect to a thermal ground 
state. In order to clarify this statement, it is useful to summarise the TBA approach. 
We assume for simplicity that the spectrum of the integrable theory consists of a single 
particle A with mass m and an exact ^-matrix S{6). In the following we consider the 
case >S'(0) = —1, which gives rise to the fermionic TBA equations. We define 

a{e) = -t\ogS{e), m = -t-^\ogS{e) . (2.1) 

The partition function at a finite temperature T and on a volume L (for L oo) is 
determined by means of the Thermodynamical Bethe Ansatz equations as follows In 
a box of large volume L, < x < L, with periodic boundary conditions, the quantization 
condition of the momenta is given by e*'^*^^'^^ Hj^i S{9i — 9j) = 1, i.e. 

niL shOi + ^ a{9i - Oj) = 27™^ , (2.2) 
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where rii are integers. Introducing a density of occupied states per unit volume pi{0) as 
well as a density of levels p{0), in the thermodynamic limit eq. ( p.2|) becomes 



2np = e + 2Tt(p * pi , (2.3) 

where e = mcosh^ and {f*g){0) = J^^d9' f {9 — 9') g {9')/ 2n. Defining the pseudo-energy 

e{9) as 

Pi 1 



(2.4) 

p 1 + e^ ' ^ ' 

the minimization of the free-energy with respect to the densities of states leads to the 
integral equation 

e = eR-(p*log{l + e-') , (2.5) 
and the partition function is then given by 

d9 



Z{L, R) = exp 



/diU 
— ch6' log (1 + e" 
27r ^ 



■e{e) 



(2.6) 



As shown in the interesting point is now that the above partition function can be 
interpreted as the one of a gas of fermionic particles but with energy given by e{9)/R. 
Namely, there is a one-to-one correspondance between the above expression ( p.6|) and 
the partition function computed according to the following thermal sum 

Z(L,R)-^^y^'-'^iO.-0.\e.-0.)Ue-''^> , (2,7) 

where the scalar products of the states are computed by applying the standard free 
fermionic rules. The above equality implies that all physical properties of the system can 
be extracted by employing the quasi-particle excitations above the TBA thermal ground 
state. Since this differs from the usual (zero-temperature) ground state, it is not sur- 
prising that its excitations do not satisfy the standard dispersion relations e = m cosh j3, 
p = msinh/?, rather they have dressed energy e = e{9)/R and dressed momentum k{9): 

e{9) = e{9)/R, k{9) = k{9) + 27f {a * pi) (9) . (2.8) 

In this contest, the rapidity 9 plays the role of a variable which simply parameterises the 
dispersion relation of the quasi-particle excitations and their S'-matrix, which is assumed 
to coincide with the original S{9i — 9j). 

The TBA allows us to compute the finite temperature one-point function of the trace 
of the stress-energy tensor [Q. In fact, we have 

TO«-(r;)o = f ^[i?i^^(i?)] , (2.9) 
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where E{R) = — log Z/L. This can be also expressed as 

{Tpn - (T;)o = mj dO-^^^ (^dne did - ^dee sh^) , (2.10) 

where the functions due and dee satisfy hnear integral equations which can be easily 
solved. The final result reads 



where 



\n=l • 



n ^fme-<'^^ 
i=i 
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)'--^(9„-i,„)ch(ei„) , (2.11) 



Let us consider now the calculation of the finite temperature one-point functions (the 
only ones which we consider in this paper). According to LeClair and Mussardo, this 
correlator is given by 

oo ^ ^ 



(f(-.')>H = i:;^7^/ 



.1=1 



(^„---^^l|O(0)|^^l---^„)eonn, (2.13) 



where the connected Form Factor of the operator O is defined as 

{On--- ei\o\e[ ■ ■ ■ e'j,,^,, = fp Qim (o|o|^„ + ^tt + ^t^^, ...,ei + t7i + t7]i,ei,...,en) 

(2.^14) 

FP in front of the above expression means taking its finite part, i.e. terms proportional 
to [l/rjiy, where p is some positive power, and also terms proportional to fji/rij, i j 
are discarded in taking the limit. With this prescription the resulting expression is an 
universal quantity, i.e. independent of the way in which the above limits are taken. 

It is easy to see that within this formalism, the finite temperature one-point function 
of the trace of the stress-energy tensor exactly coincides with its expression provided by 
the TBA, eq. (|2.11|) . In fact, the connected matrix elements of this operators are given by 

(^|T;|^)conn = 27rm2 ; (2.15) 
and by an inductive application of the form factor residue equations 

(^„---ei|T;|^^l---0„)eonn = 2777712 0(^12)0(^23) , " " " 0(^„-l,n) ch(M (2-16) 

+ permutations 

where 6'^ = 9i — 9j. Once inserted into eq. (|2.13|) , the above series coincides with the one 
of eq. ( CT) . 

In conclusion, the formalism by LeClair and Mussardo predicts, at least for the partic- 
ular thermal one-point function of T^, an exact matching with the expression determined 
by the TBA. 
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3 Delfino's Formalism 



This formalism, discussed in 0, only employs the Form Factor Approach. The finite 
temperature effects are taken into account by defining the theory on a cylinder infinitely 
extended in the space direction and a width i? = 1/T in the other direction. The particles 
entering the thermal sum are the asymptotic states satisfying the standard dispersion 
relations e = m cosh jS, p = m sinh f3 and the contribution of the n-particle asymptotic 
state to Tr[Oe^^^] is given by 

d^{R) = Jde,... dOnF^^iOn, . . . , ^il^i, . . . , 0„)e-^"^ (3.1) 

with En = m J2i=i cosh 9i and 

Define 

oo 

A^) = E^n(i?) , (3.2) 

n=0 

and normalise the thermal sum with respect the identity operator / 

^^^^ - ■ ^^-^^ 

In his paper 0] , Delfino considered for the finite part of the Form Factors entering eq. (|3.1| ) 
the symmetric limit 

TZ{dl, ...,en) = lim Fo^2n(^l +tTT + lV,...,e,+lTT + l7],e,,...,9n) , (3.4) 

and he also showed that the singular disconned parts of the Form Factors of the local 
operator O only enter through the constant factor 5'(0). All other singular terms cancel 
in the ratio ( |3.3|) . Finally, he proposed for the finite temperature one-point function the 
expression 



°° 1 1 f 

where 



Y[de,giO^,R)e 



mR cosh t 
.1=1 



TZ{0i,---.0n) , (3.5) 



g{e,R) - ^ _ ^^Q^^_„^,o,he • (3-6) 
The above formula has to be contrasted with the one given by eq. ( p.l3|) . 
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4 Main Differences and Open Problems 



There are two main differences between the two formahsms: 

• LeClair-Mussardo formahsm employs the quasi-particle excitations with respect to 
the thermal vacuum and therefore the pseudo-energy e{6), solution of the integral 
equation ( |2.3| ), whereas Delfino's formalism employs the standard asymptotic par- 
ticles at zero-temperature with energy e = m cosh 9 and momentum p = m sinh 6. 
These different choices of excitations seem somehow related to the boundary condi- 
tions adopted by the two formalisms along the space direction, i.e. in the LeClair- 
Mussardo approach one considers a box of large volume L, with periodic b.c, in the 
limit L oo, whereas in the Delfino approach one directly considers the infinitely 
extended line. Notice, however, that there is no dependence on L in the final ex- 
pressions ( |2.13| ) and ( |3.5| ) and therefore it is not a-priori clear the role played by 
the boundary conditions in thermal effects and which of the two is the appropriate 
one. 



the Form Factors entering equation ( 2.13| ) are computed according to the prescrip- 



tion given by eq. ( p.l4| ) whereas those entering equation (p.5|) are computed accord- 
ing to the symmetric limit ( |3.4D . The two different prescriptions for the finite part 
of the Form Factors produce, of course, two different results. In the case of the 
trace of the stress-energy tensor, for instance, there is already a difference for the 
two-particle Form Factor entering the thermal sum: by using the symmetric limit, 
in fact we have 

{02, 9^\Tl^\9^, 02) = 87lm'(Pi9^ -02) cosh' , (4.1) 
to be contrasted with eq. ( p.l5| ), obtained by using the other prescription. 
It is therefore evident that the two formulas, eq. ( p.l3| ) and eq. (|3.5| ), proposed for the 



one-point function at finite temperature, deeply differ in their physical justifications and 
in their technical details. To judge which of the two is the correct one it seems necessary 
to reach a better understanding of the physical principles ruling the thermal effects in 
quantum field theories. Given the present ignorance about these principles, it is therefore 
difficult to decide a-priori in favour of one or the other of the two formulations and the 
best thing one can do is to to perform some checks. Those already done and discussed in 
the literature are unfortunately inconclusive. Lukyanov @], for instance, computed the 
thermal one-point functions of the vertex operators in the Sinh-Gordon model by per- 
forming the path integral of the model and he showed that these quantities coincide with 
the ones computed in the formalism by LeClair-Mussardo. Unfortunately, the perturba- 
tive order at which he performed the computation does not permit to decide about their 
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general validity. On the other hand, Delfino |^ showed that his formalism is able to repro- 
duce the one-point function of up to the two-particle contributions but unfortunately 
he did not prove the complete equivalence of his formula with the TBA expression. 

Given the present unsatisfactory status about the validity of the two formalisms it 
is highly desiderable to perform additional checks, in particular by comparing their pre- 
dictions against a quantity which can be determined by an independent method. These 
considerations naturally select the one-point function of the trace of the stress-energy 
tensor as a check quantity for the two formulas, since its expression (|2.11|) is indepen- 



dently determined by the TBA. Hence, we have to see whether or not Delfino's formula 
reproduces the TBA result, not only up to the two-particle contribution, but also to 
higher orders (as shown above, the formula by LeClair-Mussardo coincides with the for- 
mula of the TBA). We have then two possibilities: (i) the formula proposed by Delfino is 
unable to reproduce the TBA result at higher orders; (ii) the formula proposed by Delfino 
reproduces the TBA result, alias it is just a different organization of the terms entering 
both the thermal sum and the integral equations of the TBA. In the first case, the failure 
of this check is already enough to decide about the general validity of the thermal expres- 
sions proposed by Delfino. In the second case, there would be still open the problem which 
of the two formalisms is the correct one, since their coincidence for the particular case of 
the stress-energy tensor is not expected to occur for other operators. Luckly enough, it 
is the first possibility that happens. To show the discrepancy of Delfino's formula with 
the TBA, we compare the thermal expression of the stress-energy tensor of a particular 
model which can be analytically solved. 

5 A Simplified Model 

The main technical difficulty for comparing Delfino's expression of (T^)k with the analo- 
gous expression coming from the TBA consists in solving the integral equation ( p. 31) . We 
can simplify this step by taking a local kernel, i.e. we consider an integrable model for 
which 

- 02) = 2n5{9^ - 62) . (5.1) 



For the associate S'-matrix we have 

Sm^l' (5.2) 
^ ^ I -1 if e = o . 

The integrable model defined in this way may be regarded as the limit g ^ of the Sinh- 
Gordon model. In fact, with the notation of ref. [D], the S-matrix of the Sinh-Gordon 
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model is given by 



sinh 9 — i sin 



^B(g) 



Ssh{0) = "-^^ J, . , (5.3) 

sinh 6 + i sin ^ 

with B{g) = g^^g2 ■ It is convenient to define B{g) = 2a. For the corresponding kernel 
we have 

2 sin na cosh 9 
sinh^ 9 + sin^ na ' 

and in the limit a ^ we have 



shiO) = . ,2. , .2 > (5-4) 



Iimjshi0) = 2n5i9) . (5.5) 



By using the kernel ( |5.1j ), the integral equation ( |2.3| ) becomes 

£(^) = mR cosh ^ - ln(l + e''^'^^) , (5.6) 

whose solution is given by 

6(9) = In (^e™^-°«i^« - i) . (5.7) 

Hence 

f(9) e-^(^) = ' = e— , (5.8) 
1 + 

and inserting into the TBA formula ( |2.11| ), we have 

(Wfi - (r^)0 = 27rm2 / L-mfi cosh e^g-2miJ cosh e^g-3mi? cosh e^... ^ (gg) 

J-oo 27r L J 



For the purpose of comparing with Delfino's prediction, it is convenient to leave explicitly 
the n-particle contributions to the thermal average, although it is evident that the above 
series can summed to 

{Tll)n - (T;)o = 2.m- , (5.10) 

which is nothing else but the thermal one-point function of for a free bosonic theory. 

Let us consider now the Form Factors of associated to the simplified model with 
kernel (|5.1|) . In virtue of the observed equivalence of this theory with a particular limit of 
the Sinh-Gordon model, the Form Factors can obtained by a careful (7^0 limit of the 
corresponding quantities of the Sinh-Gordon model. They were computed in and can 
be expressed as 

(0|T;(0)|^,,...,^„) = -^^-^U^^ Q„(a;i,...,x„)nf^ • (5.11) 

Few words on the above expression. The explicit form of Fmin{9) can be found in ||^. For 
our purposes we only need the functional equation satisfied by Fmin{9) 

sinh 9 

F^in{9)Fmin{9 + m) = . (5.12) 

smh -|- « sm vra 
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Qn is a symmetric polynomial in the variables Xj = e^' given by 

Qn{xi, ...,Xn) = detMij , 
with the (n — 1) x (n — 1) matrix Mij given by 

Mij = 02i^j[i - J + 1] . 

In the above equation the symbol [n] is defined by 

sin(na) 

n = — , 

sm a 

and (Tfc is the elementary symmetric polynomial given by the generating function 

n n 

Y[{x + Xi) = x""''(Tfc(xi, a;2, ...,Xn) . 

i=l k=0 

In the limit a — 0, the first polynomials Q„ are given by 
Q2 = cTi ; 

Q4 = 0-10-20-3 ; (5.13) 
Qe = cTiCTs [f^2cr30-4 + Sorgo-Q - 4:{aia2crQ + 0-40-5)] . 



5.1 Two— particle contribution 



By using eq. (|5.11|), let us compute 



{e2,ei\TI::\eu02) = Hm limlOmei+t-K + r],, 62 + t7c + r]2, 61,62) . (5.14) 

We will consider the contributions coming from the different terms in ( p. Ill ) separately. 

By using the functional equation (|5.12|) , for the product of Fmin{6ij) we have, in the 
above limit ^ 

j<-j smn ui2 + sm vra 

For the polynomial of the denominator we have 



]^(xi + a;j) — ^ AiA2XiX2{xi + X2)'^{xi - X2y , (5.16) 

i<j 

where = (1 — e*'"') ~ — ir^fc. Finally, for the polynomial Q4 in the numerator we obtain 

Q4 — ^ XiX2{xj + xl)[{Al + Al)xiX2 + AiA2ixl + xD] . (5.17) 

We have now two possibilities. The first consists of keeping in the above expression 
only the term multiplying the combination ^1^12 (and disregarding those multiplying 
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{Af + A2)). This leads to the computation of the connected Form Factor. In this case, 



combining all terms and taking the limit (|5.14|) , we have 



2 / 2 sin 7ia cosh 612 



h,Oimei, e2)conn = 47rm^ —-2^ cosh ^12 . (5.18) 

V smh 612 + sm vra / 



By taking now the limit a ^ and using eq. (p.5|), we have 

(^2,^l|7;1^1,^2)conn = 471771^(^1 " ^2) COsh ^12 , (5.19) 

in agreement with eq. ( p.l5|) . 

The second possibility consists of taking the symmetric limit considered by Delfino. 
This is obtained by taking Ai = A2. In this case, the symmetric limit of eq. ( |5.14| ) 
produces 

/a a q\ q 2 sin vra cosh ^12 \ 2^12 

{d2,di\T^\di,02)sym = S7Tm — — ^osh — . (5.20) 

^ \smh 6'i2 + sm vray 2 

By taking now the limit a ^ and using eq. (^.5|), we have 

(^2,^i|T;|^i,^2)sym = 87rmV(^i - ^2) cosh^ ^ . (5.21) 



Let us consider the expression ( |3.5| ) up to the the two-particle contribution. For the 
function g{6, R) we have 

fi'(^!-R) — ^ _ ^^g^g-_m_Rcosh6» ~ ^ + e~^^^°^^^ ' (5.22) 



and then 

/+OO rjf) 
-00 Z7T 



g— mi? cosh 9 g— 2m_Rcosh( 



(5.23) 



+ 2 h ■ ■ 

]^ _j_ g— m_Rcosh(? _|_ g— »n_Rcosh6^2 

Expanding this expression in power of e"™-^™'^'^^ up to e"^™-^™'^'^^ we have 

/+OO rlQ y -, 
_ g-miJ cosh e^g-2m/? cosh e^(>^^g-3m/? cosh _ ^5 24) 

-00 27r '- 

Comparing now this expression with eq. ( [5.8| ), we explicitly confirm the agreement found 
at this order by Delfino in his paper. 

5.2 Three— particle contribution 

By using eq. (|5.11| ), let us compute 



93,^2,^i|T;|^i,^2,^2) = liin lim lim (0|T;|^i+m+7/i, ^s+^vr+r^s, ^3+^vr+7/3, ^i, ^2, ^^3) • 

(5.25) 
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As before, let us consider the contributions coming from the different terms separately. 
By using the functional equation (|5.12|) , for the product of Fmin{Gij) we have, in the above 
limit 



sinh^ 612 \ ( sinh^ 6'i3 \ / sinh^ ^23 



sinh 6'i2 + sin vra / \sinh 6'i3 + sin vra/ \sinh 6'23 + sin Txa^ 

(5.26) 

For the polynomial of the denominator we have 

+ — > ^1^2^3X1X2X3 (a^i — a;2)(a;i — X3)(x2 — X3) = 

X1X2X3 

)^(sinh^i2sinh^i3sinh^23)^ • (5.27) 

For the polynomial Qq, we have two possibilities. The first consists of keeping only the 
term multiplying the combination A1A2A3 (and disregarding all other expressions which 
multiply the other monomials like Af, AIA2 etc.). This leads to the computation of the 
connected Form Factor. In this case we have 

gconn ^ AiA2A^XiX2X^{xl + xl){xl + xDixj + xj) X (5.28) 

(xiX2)^(Xi + X2 - 2X3) + (xiX3)^(Xi + X3 - 2X2) + (X2X3)^(X2 + X3 - 2Xi) 

and for the connected Form Factor, combining all terms, we obtain 

/ft ft ft iT^^lft ft ft\ - 0-rrr,2/^ 2 siu TTQ COsh ^12 \( 2 siu TTg COSh ^23 \ 
(6^3,fc^2,t^l|J^TI^l>^2,t^2)conn - 27rm . 2n \ ~ ■ 1 2 n ] — P< 

\ smh t7i2 + sm ira J \ smh (723 + sm ira J 
sinh^ ^13 
sinh^ 6*13 + sin^ na 

By taking the limit a — > 0, we obtain the result reported in formula ( p.l6|) . 

The second possibilities consists of considering the symmetric limit, which is obtained 
by taking A = Ai = A2 = A^. Other terms enter the polynomial in this case and we have 

(xiX2)^(x? + xl- 2x1) + {xix^fixj + xl- 2x1) + {x2X3)^{xl + xl- 2xj) + 

-2X1X2X3 [X2(xi - X3)^ + Xi(x2 - X3)^ + X3(xi - X2)^] ] (5.30) 

Therefore, combining all the different contributions, in the symmetric limit we have 

/n n n iT^ft ft ft\ - O-rrr,^ 2 siu TTQ COSh ^12 \( 2 siu TTQ COsh ^13 \ 
(6^3,6^2,t^l|j;rFl5^2,t^3)symm - 2TTm . 2n \ ~ ■ ,2n ] — P< 

\smh t7i2 + sm rca J \smh t^i3 + sm rra / 
sinh^ 6 

^ . , 2 n ^T-o [2 (cosh 6^12 + cosh 6*13 + cosh6'23) + 3] x 

smh 6^23 + sm na 

cosh ^ cosh %i (2 cosh^ ^ - 1) 
^ n 1 n ^ permutations . 



X 2 ; — ~ cosh 6'i3 + permutations . (5.29) 



cosh 6'i2 cosh 6*13 cosh % 
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By taking now the limit a ^ and using eq. ( p.5| ), we obtain 

CO 

cosh 9i2 cosh 6^13 



cosh — cosh ^ 

^3,e2,e,me,,e2,e2)sym = 27rm20(ei - ^2)^ - ^3) — ttt — x 



(2 cosh^ ^ — 1) 

[2(cOsh^l2 + cosh + COSh^^23) + 3]^ ^^^^^^ ' 

+permutations . (5.31) 



Once inserted into eq. ( |3.5| ), we have 

(5.32) 

and by making an expansion up to e"^™-^™*^"^^ we have 



-oo 2ti 



g— mi? cosh 6 g— 2m_Rcosh9 g g— SmRcoshS 

h 2 1 h ■ 

_j_ g— mRcoshS _j_ g— mi?cosh6^2 2 (1 -|- g^^-^coshS^S 



/-+00 r 

(^/T/fl - l^/TJo = 27rm TT ^ +^ + + C;(e j 

J —oo ZTT L Z 

(5. 

i.e. the third order coefficient disagrees with the corresponding coefficient of eq. (p.9|). 



6 Conclusions 

In this paper we have critically analysed the status of the thermal formalism for two- 
dimensional integrable field theory by comparing the approach proposed by LeClair and 
Mussardo with the approach proposed by Delfino. Whereas the ffist approach is able 
to reproduce the one-point function of as given by the TBA, the second one is in 
agreement with the TBA formula only up to the two-particle contribution and differs 
otherwise. This has been explicitly shown by considering a simple integrable model, 
where all calculations can be performed analytically without relying on the solution of 
integral equation. It would be useful to further explore the subject and see whether or 
not the approach by LeClair and Mussardo passes other tests. 
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